Microscopic sub-barrier fusion calculations for the neutron star crust 
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Fusion of very neutron rich nuclei may be important to determine the composition and heating 
of the crust of accreting neutron stars. Fusion cross sections are calculated using time-dependent 
Hartree-Fock theory coupled with density-constrained Hartree-Fock calculations to deduce an ef- 
fective potential. Systems studied include I'^O-F^'^O, ^^O+^^O, ^^O+'^'^O, ^"^C+^^O, and i^C-h^^O. 
We find remarkable agreement with experimental cross sections for the fusion of stable nuclei. Our 
simulations use the SLy4 Skyrme force that has been previously fit to the properties of stable nuclei, 
and no parameters have been fit to fusion data. We compare our results to the simple Sao Paulo 
static barrier penetration model. For the asymmetric systems ^^C-|-^*0 or ^''0+'^*0 we predict 
an order of magnitude larger cross section than those predicted by the Sao Paulo model. This is 
likely due to the transfer of neutrons from the very neutron rich nucleus to the stable nucleus and 
dynamical rearrangements of the nuclear densities during the collision process. These effects are 
not included in potential models. This enhancement of fusion cross sections, for very neutron rich 
nuclei, can be tested in the laboratory with radioactive beams. 



PACS numbers: 25.60.Pj, 26.60.-I-GJ, 97.80.Jp 

I. INTRODUCTION 

There is broad interest in the sub-barrier fusion of neu- 
tron rich nuclei, in astrophysics as we discuss in Sec. |I A[ 
in the laboratory as discussed in Sec. |IB[ and in theory, 
see Sec. ITCl 



A. Fusion in accreting neutron stars 

Neutron stars, in binary systems, can accrete mate- 
rial from their companions. This material undergoes a 
variety of nuclear reactions. First at low densities, con- 
ventional thermonuclear fusion takes place, see for ex- 
ample mil]. Next at higher densities, the rising electron 
Fermi energy induces a series of electron captures [3^ to 
produce increasingly neutron rich nuclei. Finally at high 
densities, these very neutron rich nuclei can fuse via py- 
cnonuclear reactions. Pycnonuclear fusion is induced by 
quantum zero point motion [4j . The energy released, and 
the densities at which these reactions occur, are impor- 
tant for deternrining the temperature and composition 
profile of accreting neutron star crusts. 

Superbursts are very energetic X-ray bursts from ac- 
creting neutron stars that are thought to involve the un- 
stable thermonuclear burning of carbon [316]. However, 
the best current simulations do not reproduce the con- 
ditions needed for carbon ignition because they have too 



low temperatures [7]. An additional heat source, from 
fusion or other reactions, could raise the temperature 
and allow carbon ignition at densities that reproduce ob- 
served burst frequencies. Alternatively, there could be a 
very low energy resonance in the ^^C -^^C fusion cross- 
section that could explain superburst ignition [S]. 

Haensel and Zdunik have calculated pycnonuclear fu- 
sion reactions at great densities in the inner crust of neu- 
tron stars [9] using a very simple crust composition con- 
sisting of a single average nucleus. See also the work by 
Satolin]. Gupta et al. mention the potential impact of 
electron capture and neutron emission cascades in their 
model on fusion |TT] . Here a cascade of multiple electron 
captures and neutron enrissions can reduce the atomic 
number of nuclei until the lower Coulomb barrier allows 
rapid pycnonuclear fusion. 

Horowitz et al. |12| calculate the enhancement in fusion 
rates from strong ion screening using molecular dynam- 
ics simulations, and find that ^'^O -t- ^*0 can fuse near 
10^^ g/cm'^, just before neutron drip. Lau uses a reac- 
tion network to follow the composition of an accreting 
fluid element in the neutron star crust |13] . She con- 
siders a range of fusion reactions involving neutron rich 
isotopes of chemical elements that include carbon to mag- 
nesium [13]. These reactions tend to occur at densities 
just above neutron drip (> few xlO^^ g/cm"^) and be- 
yond. Therefore we expect fusion reactions of neutron 
rich isotopes, near the drip line, to be important. Fur- 
thermore, this fusion can take place in the background 
neutron gas that is present in the inner crust of a neutron 
star. 
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B. Fusion experiments with neutron rich 
radioactive beams 



Recent advances in radioactive beam technologies have 
opened up new experimental possibilities to study fusion 
of neutron rich nuclei. Furthermore, near barrier fusion 
cross sections are relatively large so experiments are fea- 
sible with modest beam intensities. For example, at the 
GANIL-SPIRAL facility a reaccelerated beam of ^°0 was 
used to measure near- and sub-barrier fusion of ^'^O on 
""^^C [15]. In addition, measurements are possible at the 
TRIUMF ISAC facility and in the near future at the 
NSCL ReA3-6 reaccelerated beam facility. Note that the 
dynamics of the neutron rich skin of these nuclei can en- 
hance the cross-section over that predicted by a simple 
static barrier penetration model. For example, neutrons 
may be transferred from the neutron rich beam to the 
stable target. Fusion of very neutron rich nuclei, near 
the drip line, raise very interesting nuclear structure and 
nuclear dynamics questions. 



C. Calculations of sub-barrier fusion cross sections 

In the absence of a practical ab-initio quantal many- 
body theory for sub-barrier fusion all approaches involve 
the calculation of an ion-ion potential barrier, usually as 
a function of the nuclear separation coordinate R, and 
the solution of the corresponding one-body Schrodinger 
equation for the transmission probability and the fu- 
sion cross-sections. Many of the phenomenological and 
semi-microscopic potentials for fusion utilize the double- 
folding method (TBI [TTj , which is based on the physical 
assumption of frozen densities or the sudden approxima- 
tion. As the name suggests, in this approximation the nu- 
clear densities are unchanged during the computation of 
the ion-ion potential as a function of the internuclear dis- 
tance. While asymptotically fusion potentials may be de- 
termined from Coulomb and centrifugal interactions, the 
short distance behavior strongly depends on the nuclear 
surface properties and the readjustments of the combined 
nuclear system, resulting in potential pockets, which de- 
termine the characteristics of the compound nuclear sys- 
tem. For this reason approaches based on frozen densities 
start to fail as the nuclear surface overlap begins to be 
substantial due to lack of rearrangements and exchange 
effects. This is the situation for fusion at deep sub-barrier 
energies for which the inner turning point of the poten- 
tial barrier happens at a large density overlap. In the 
coupled-channel approach this has been addressed with 
the addition of a repulsive core potential at small nuclear 
separations [HllTlj. 

During the past several years, we have developed 
a new microscopic approach for calculating heavy-ion 
interaction potentials that incorporates all of the dy- 
namical entrance channel effects included in the time- 
dependent Hartree-Fock (TDHF) description of the col- 
lision process j20j. The method is based on the TDHF 



evolution of the nuclear system coupled with density- 
constrained Hartree-Fock calculations (DC-TDHF) to 
obtain the ion-ion interaction potential. The formalism 
was applied to study fusion cross-sections for the sys- 
tems 132Sn-H64]v^i ^ ^^m+^^m 122j, 16O+208pb [^, 

I32,i24gj-^^962j. [2^^ a^jjjj ^;o reactions leading to the super- 
heavy element Z = 112 [53], among others. In all cases 
a good agreement between the measured fusion cross- 
sections and the DC-TDHF results were found. This is 
rather remarkable given the fact that the only input in 
DC-TDHF is the nuclear effective interaction, and there 
are no adjustable parameters. 

There is a great deal of experimental information on 
low energy fusion cross-sections for light stable nuclei 
such as i^C gillTI and ^^o ^,^. For these nuclei, 
barrier penetration models work well for energies near 
the Coulomb barrier ^30J. However, recently Jiang et al. 
discussed fusion hindrance at extreme sub coulomb bar- 
rier energies ^3lj- Much less information is available for 
the fusion of very neutron rich light nuclei. 

The Sao Paulo model of fusion calculates an effective 
potential based on the density overlap between colliding 
nuclei 32J. Sub-barrier fusion cross-sections can then be 
calculated via tunneling. The model can be easily applied 
to a very large range of fusion reactions and qualitatively 
reproduces many experimental cross-sections [33| . Re- 
cently this model was used to tabulate astrophysical S 
factors describing fusion of many carbon, oxygen, neon 
and magnesium isotopes for use in astrophysical simula- 
tions [34 , see also Ref. |35j . 

This paper is organized as follows. Our time depen- 
dent Hartree Fock formalism is discussed in Secjllj while 
Sec. |III| presents results for effective potentials and fusion 
cross sections. Finally, these results are discussed and we 
conclude in Sec. IIVI 



II. FORMALISM 

We describe our time-dependent Hartree-Fock based 
formalism in Sec. |II A| and then, for comparison, re- 
view the simple Sao Paulo barrier penetration model in 
Secimi 



A. TDHF and DC-TDHF 

In this paper, we utilize the DC-TDHF method for 
the calculation of the ion-ion potential barriers. It is 
generally acknowledged that the TDHF method provides 
a useful foundation for a fully microscopic many-body 
theory of low-energy heavy-ion reactions ^36] based on a 
time-dependent nuclear energy-density functional. Dur- 
ing 1970's and 1980's the TDHF theory has been widely 
used in the study of fusion excitation functions, deep- 
inelastic scattering of heavy mass systems, and nuclear 
molecular resonances, while providing a natural founda- 
tion for many other studies. An account of some of the 



earlier TDHF applications can be found in Refs. '3S'j '37] . 
With modern supercomputers it has become feasible to 
carry out very accurate nuclear structure and reaction 
studies in full three-dimensional space and using the 
complete form of the nuclear energy-density functional. 
The DC-TDHF approach is analogous to calculating mi- 
croscopic potential energy surfaces with the constrained 
Hartree-Fock method using e.g. quadrupole and octupole 
constraints. However, in this approach the TDHF time- 
evolution takes place with no restrictions. At certain 
times during the time evolution we perform a paral- 
lel static Hartree-Fock minimization while holding the 
neutron and proton densities constrained to be the cor- 
responding instantaneous TDHF densities. In essence, 
this gives us the TDHF dynamical path in relation to 
the multi-dimensional static energy surface of the com- 
bined nuclear system. The DC-TDHF approach [20| pro- 
vides the means for extracting ion-ion potentials from the 
TDHF evolution of the nuclear collision as follows: The 
ion-ion interaction potential is given by 



Vuc{R) = Eucifi) - Ea, - Ea, 



(1) 



where -Edc is the density-constrained energy at the in- 
stantaneous separation between nuclear centers R{t), 
while -Eai and £'a2 ^re the binding energies of the two 
nuclei obtained with the same energy-density functional. 
The TDHF evolution also provides us with a coordinate- 
dependent mass, M{R), which should be used in solving 
the resulting Schrodinger equation, using the conserva- 
tion of energy 



M{R) 



2[E, 



vocm 



i?2 



(2) 



where the collective velocity R is directly obtained from 
the TDHF evolution. The ^-dependence of this mass at 
lower energies is very similar to the one found in con- 
strained Hartree-Fock calculations '38' with a constraint 
on the quadrupole moment. On the other hand, at higher 
energies the coordinate dependent mass essentially be- 
comes flat, which is again a sign that most dynamical ef- 
fects are contained at lower energies. The peak at small 
R values is due to the fact that the center-of-mass en- 
ergy is above the barrier and the denominator of Eq. ([2]) 
becomes small due to the slowdown of the ions. 

All of the dynamical features included in TDHF are 
naturally included in the DC-TDHF potentials. These ef- 
fects include neck formation, mass exchange, internal ex- 
citations, deformation effects to all order, among others. 
Alternatively, instead of solving the Schrodinger equa- 
tion with coordinate dependent mass parameter M{R) 
for the heavy- ion potential V-aciR), we can instead use 
the constant reduced mass /i and transfer the coordinate- 
dependence of the mass to a scaled potential V{R) using 
the well known exact point transformation [231 138) 



di?=(ffl) dR 



(3) 



The potential V{R), which includes the coordinate- 
dependent mass effects differs from the T^dc(-R) only in 
the interior region of the barrier. Further details can be 
found in Ref. E31. 
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FIG. 1: (Color online) Snapshots of the nuclear density con- 
tours in the x—z plane during the TDHF time-evolution of the 
^^0-|-^^0 system at a collision energy -Etdhf ~ 12 MeV are 
shown for four different ion-ion separation values of R. The 
top panel corresponds to the peak of the potential barrier at 
10.0 MeV, the subsequent panels show the nuclear density at 
the inner turning point of the potential barrier corresponding 
to the cm. energies, 9.4, 6.5, and 3.7 MeV, respectively. 

The fusion barrier penetrabilities Ti(£'c.m.) are ob- 
tained by numerical integration of the two-body 
Schrodinger equation 



-h"^ (P LiL + l)h^ 



2^ di?2 



2^i?2 



V{R)-E 



^ = , (4) 



using the incoming wave boundary condition (IWBC) 
method [39]. IWBC assumes that once the minimum 
of the potential is reached fusion will occur. In prac- 
tice, the Schrodinger equation is integrated from the po- 
tential minimum, i?min, where only an incoming wave 



is assumed, to a large asymptotic distance, where it 
is matched to incoming and outgoing Coulomb wave- 
functions. The barrier penetration factor, T^iyEcm) is 
the ratio of the incoming flux at i?min to the incoming 
Coulomb flux at large distance. Here, we implement the 
IWBC method exactly as it is formulated for the coupled- 
channel code CCFULL described in Ref. [JD]- This gives 
us a consistent way for calculating cross-sections at above 
and below the barrier via 






(5) 



L=0 



starts with the double folding potential Vf{R) |32j . 
Vf{R) ^ I d\^d\2Pi{.ri)p2{r2)V^5{vi - rs - R) . (6) 

Here pi and p2 are the ground-state densities of the 
two nuclei and Vq = —456 MeV-fm^. From Vp a 
nonlocal potential V{R,E) is constructed, V{R,E) = 
VF{R)e~^^ 1'^ , where v is the local relative velocity f35] 
and c is the speed of light. Finally, tunneling through the 
Coulomb plus V{R-, E) potentials is calculated in a WKB 
approximation to infer the fusion cross section [4i ,32; . 




FIG. 2; (Color online) Snapshots of the nuclear density con- 
tours in the x—z plane during the TDHF time-evolution of the 
16q_|_24q gyg^gjjj g^^ g^ collislon Bncrgy i?TDHF = 12 MeV are 
shown for four different ion-ion separation values of R. The 
top panel corresponds to the peak of the potential barrier at 
8.5 MeV, the subsequent panels show the nuclear density at 
the inner turning point of the potential barrier corresponding 
to the cm. energies, 6.0, 4.0, and 3.0 MeV, respectively. 



B. Simple Sao Paulo barrier penetration model 

In this section we describe the simple Sao Paulo barrier 
penetration model to calculate fusion cross-sections. This 



III. RESULTS 

In this section we present results for effective poten- 
tials and fusion cross-sections using the two approaches 
outlined in the preceding section. 

Time-dependent Hartree-Fock calculations were done 
in 3-D geometry and using the full Skyrme interaction 
including all of the time-odd terms in the mean-field 
Hamiltonian [41] . The Skyrme parametrization used was 
SLy4 [32] ■ In addition to providing a good description 
of nuclei, this interaction has been used to produce an 
equation of state for neutron stars [43] . For the reactions 
studied here, the lattice spans 44 fm along the collision 
axis and 22 fm in the other two directions. Derivative op- 
erators on the lattice are represented by the Basis-Spline 
collocation method. One of the major advantages of this 
method is that we may use a relatively large grid spacing 
of 1.0 fm and nevertheless achieve high numerical accu- 
racy. The initial separation of the two nuclei was 18 fm 
for central collisions and the TDHF collision energy was 
-£'TDHF = 12 MeV. The time-propagation is carried out 
using a Taylor series expansion (up to orders 10 — 12) 
of the unitary mean-field propagator, with a time-step 
At — 0.4 fm/c. The accuracy of the density constraint 
calculations is commensurate with the accuracy of the 
static calculations. 

Calculations were done on a 16-processor Linux work- 
station using our OpenMP TDHF code with almost a 
100% parallel efficiency. We have performed density con- 
straint calculations every 10 — 20 time steps. For the light 
systems considered here, one full calculation takes about 
12 hours. 

In Fig. [T]we plot snapshots of the nuclear density con- 
tours in the x — z plane during the TDHF time-evolution 
of the ^^0-1-^^0 system at four different ion-ion separa- 
tion values of R. The top panel corresponds to the peak 
of the potential barrier at 10.0 MeV, the subsequent pan- 
els show the nuclear density at the inner turning point 
of the potential barrier corresponding to cm. energies of 
9.4, 6.5, and 3.7 MeV, respectively. As one can observe 
from these densities for low cm. energies the nuclei have 
a substantial overlap and develop a neck which is very 
different from two overlapping spherical densities as in 
the frozen-density approach. Figure [2] shows the density 
contours for the asymmetric ^^0-1-^3 system. Again, 



the top panel represents the density contours correspond- 
ing to the peak of the corresponding potential barrier at 
8.5 MeV, the subsequent panels show the same quantity 
at barrier energies 6.0, 4.0, and 3.0 MeV, respectively. 
The difference between the two figures is interesting. For 
the neutron-rich system the barrier peak is at a larger 
R value since the nuclei come into contact sooner due to 
the extended neutron skin of the ^*0 nucleus. 

For simplicity the Sao Paulo model often assumes 
Woods-Saxon static nuclear densities with radius pa- 
rameter R — 1.31v4^/"^ — 0.84 fm and diffuseness a — 
0.58 fm [31|32]. However there are also Sao Paulo model 
results using mean field densities [H]. In Fig. [s] we 
compare these Woods-Saxon ground state densities with 
those calculated with a spherical Hartree-Fock program 
using the Skyrme SLy4 interaction. As one can see the 
surface properties calculated with the two-approaches 
agree reasonably well, whereas the central dip in the den- 
sity, characteristic to some light nuclei, is not reproduced 
by the Woods-Saxon shape. 
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FIG. 3: (Color online) Baryon density of oxygen isotopes ^®0 
(black) and ^*0 (red) versus radial distance r. The solid lines 
are HF results for the SLy4 interaction while the dashed lines 
are simple parameterized Woods-Saxon densities used in the 
Sao Paulo barrier penetration model, see Sec. |IIB[ 



In Fi g. [4| we have plotted various quantities discussed 
in Sec. II A for the ^^0-1- ^^O system. We observe how 
the coordinate-dependent mass modifies the inner part 
of the barrier when the point transformation is per- 
formed. The dashed line shows the potential Vdc{R) 
calculated directly from TDHF with density constraint. 
The coordinate-dependent mass M{R)/fi is also shown 
on the same plot. As we can see this ratio equals one 
when the two nuclei are far apart and starts to deviate 
as the nuclei slow down and finally peaks when R ap- 
proaches zero. The point transformation Eq. ([3]) results 
in the potential V{R). We stress that the transforma- 



tion given in Eq. ([3]) is carried out for numerical conve- 
nience: The original Schrodinger equation which involves 
the coordinate-dependent mass M{R) and the potential 
Vdc (R) is mathematically equivalent to the transformed 
Eq. (HI) which contains the constant reduced mass /i and 
the transformed potential V{R). However, from a nu- 
merical point of view, Eq. Q is easier to solve because 
of the simpler structure of the kinetic energy operator. 
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FIG. 4: (Color online) DC-TDHF potential barriers for the 
^''0+ ^^O system. The dashed line shows the potential 
Vdc{R) calculated directly from TDHF with density con- 
straints. The coordinate-dependent mass M{R)/jj, is also 
shown on the same plot. The point transformation Eq. (|3| re- 
sults in the potential V{R). Also shown is the point Coulomb 
potential. 



In Fig. [5] we plot the ion-ion potentials obtained using 
the density-constrained TDHF approach for various oxy- 
gen isotopes. These potentials include all the dynamical 
effects as well as the coordinate dependent mass. We see 
that the ""^^0-1- ^^O potential (black solid curve) has the 
highest barrier of about 10 MeV. As the nuclear mass 
increases the nuclei come into contact earlier, which cor- 
responds to a larger R value, and deviate from the point 
Coulomb potential (green solid curve), consequently hav- 
ing a lower barrier height. The barrier heights in decreas- 
ing order are 10.00, 9.24, 8.54, and 7.95 MeV. The rea- 
son for potentials extending to smaller R values in the 
inner part of the barrier is due to the fact that the col- 
lision energy is effectively higher for lower barriers thus 
the system can reach a more compact shape before fu- 
sion. With the SLy4 interaction the ^^0 nucleus is barely 
bound, while the experimental drip line is at ^"^O, so we 
have included this nucleus as well. 

Excitations are believed to have a significant impact 
on the outcome of the fusion reactions. The excitations 
can range from the entrance channel quantal excitations 
of the projectile and target, as in the coupled-channel 
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FIG. 5: (Color online) Potential barriers for the 0+0 system 
obtained from density constrained TDHF calculations, ^"O-f 
^®0 (black solid curve), ^^0-1- ^''O (red short-dash curve), 
^'^0+ ^■'O (blue long-dash curve), and ^^'0+ ^®0 (orange dot- 
ted curve). Also shown is the point Coulomb potential. 
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FIG. 6: (Color online) Neutron transfer (blue) and proton 
transfer (red) to ^^O from the neutron-rich reaction partner 
versus separation R during the fusion of ^^0 on ^'*0 (solid 
lines) and on ^*0 (dashed lines). 



approach, to collective excitations of pre-equilibriuni sys- 
tem, to compound nucleus excitations. These can be fur- 
ther influenced by particle transfer, pre-equilibrium emis- 
sions, and evaporation, among others. Theoretically such 
effects are commonly introduced by hand into various re- 
action models. However, the influence of excitations on 



nuclear reaction dynamics remains to be a difficult and an 
open problem as it combines both nuclear structure and 
dynamics under nonequilibrium conditions. In Fig. [6] we 
plot the average number of neutrons and protons trans- 
ferred to the ^^O nucleus during the early stages of the 
TDHF coUisions for the i^O-f24o (solid curves) and the 
16q_|_28q ((fashed curves) systems. Neutrons are denoted 
by blue curves and protons by red curves. These trans- 
fers are calculated using the standard method in TDHF, 
which is to take a cut at the point of minimum density in 
the overlap region between the two nuclei and integrate 
the densities in the left and right halves of the collision 
box. A number of interesting things can be observed 
from this plot; first is the fact that most of the transfer 
seems to start after we pass the potential barrier peaks. 
This indicates that particle transfer primarily modifies 
the inner part of the barriers and not so much the bar- 
rier heights. The other observation is that on average 
about one neutron is transferred from the ^^O to ^^O in 
the region of R values relevant for fusion cross-sections, 
whereas two neutrons are transferred from ^^O to ^^O in 
the same region. We also note that these transfers occur 
rather rapidly. 
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FIG. 7: (Color online) Cross-section versus center-of-mass en- 
ergy for fusion of oxygen isotopes. Experimental results (filled 
circles) are from Ref. [29| 



We have obtained the fusion cross-sections for the 
0-f O systems by numerical integration of Eq. J4| to cal- 
culate the transmission probability and Eq. (Isr for the 
fusion cross sections. The resulting cross-sections are 
shown in Fig. [t) We observe that for the ^^0-1-^^0 sys- 
tem barriers we obtain a very good description of the low- 
energy experimental fusion cross-sections. The higher en- 
ergy part of the fusion cross-sections are primarily deter- 
mined by the barrier properties in the vicinity of the bar- 
rier peak. On the other hand sub-barrier cross-sections 
are influenced by what happens in the inner part of the 



barrier and here the dynamics and consequently the co- 
ordinate dependent mass becomes very important. Also 
shown in Fig.JTJare the fusion cross-sections for ^^0+^*0, 
160+280^ and the 240-^2^0 systems. The behavior of 
these cross-sections can be deduced from the potential 
barriers of Fig. [5) The highest fusion cross-sections be- 
long to the collision of ^^O with the most neutron rich 
isotope ^^0. The interesting observation is that the col- 
lisions involving ^^O and a neutron-rich isotope seem to 
have a larger fusion yield than the ^'*0+^'*0 system. 




FIG. 8: (Color online) Potential barriers for the C-l-O system 
obtained from density constrained TDHF calculations, ^^C-|- 
^''O (black solid curve) and ^^C-h ^*0 (blue dashed curve). 
Also shown is the point Coulomb potential (solid green curve). 



In Fig. [8] we show the DC-TDHF potential barriers for 
the C+0 system. The higher barrier corresponds to the 
^^G+ ^^O system and has a peak energy of 7.77 MeV. The 
barrier for the ^^C+ '^"'O system occurs at a slightly larger 
R value with a barrier peak of 6.64 MeV. Figure |9] shows 
the corresponding cross-sections for the two reactions. 
Also shown are the experimental data from Ref . [3T] . The 
DC-TDHF potential reproduces the experimental cross- 
sections quite well for the ^^C-|- ^^O system. Again the 
cross section for the neutron rich ^^C+^^0 is seen to be 
larger than that for ^^C +^^0. 

Some of the strong energy dependence of the fusion 
cross section (7{E) can be taken into account with the 
Astrophysical S factor 



S{E) ^a{E)EExp[2Trr]] 



(7) 



at center-of-mass energy E and Sommerfeld parameter 
T] = ZiZ2e^ /{hv). Finally, the relative velocity of the 
nuclei is u = \J2EI [i for a system of reduced mass /i. 

The S factor for our TDHF calculations is shown 
in Figs. 10 and 11 Also shown are S factors from 
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FIG. 9: (Color online) Cross-section versus center-of-mass en- 
ergy for fusion of ^^C with oxygen isotopes. Experimental 
data (circles) are from Ref. |31| . 



accurate nine parameter fits to Sao Paulo model re- 
sults |M| (dashed lines). In Fig. [lO] the TDHF S fac- 
tor for ^^0-1- "'^^O agrees well with data. Presumably a 



small change in the Sao Paulo potential will allow the 
Sao Paulo results to also agree well with data. How- 
ever, there are interesting differences for more neutron 
rich systems. We note that the Sao Paulo S factor for 
the symmetric system ^*0-|-^'*0 is larger than the TDHF 
result while the Sao Paul S for the asymmetric system 
16q^24q jg significantly smaller than the TDHF result. 
We believe this pattern is due to dynamics present in the 
TDHF formalism that is not included in the simple Sao 
Paulo model. 

Figure [Tl] shows a similar pattern for reactions in- 
volving ^^C. The TDHF results agree well with data for 
^^C-t-^^0 and predict a much larger S factor, compared 
to the Sao Paulo model, for the neutron rich asymmetric 
system ^^C-|-^'^0. We discuss these differences between 
our TDHF results and Sao Paulo model results in the 
next section. 

Finally, fusion in accreting neutron stars may take 
place in the inner crust where there is a background gas 
of neutrons. This background gas could impact fusion 
cross sections. To investigate this we are performing 
TDHF simulations where the initial conditions involve 
both the two nuclei, appropriately boosted towards each 
other, and some unbound neutrons, that are uniformly 
distributed in the simulation volume. We will report 
these results in a later publication. 
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FIG. 10: (Color online) Astrophysical S factor versus center 
of mass energy for fusion of oxygen isotopes. Solid lines show 
TDHF results while dashed lines are for the Sao Paulo barrier 
penetration model. Experimental results (circles) are from 
ref. ESI. 
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FIG. 11: (Color online) Astrophysical 5* factor versus center of 
mass energy for fusion of ^^C with oxygen isotopes. Solid lines 
show TDHF results while dashed lines are for the Sao Paulo 
barrier penetration model. Experimental results (circles) are 
from ref. 1311. 



IV. SUMMARY AND CONCLUSIONS 

Fusion of very neutron rich nuclei may be important 
to determine the composition and heating of the crust of 
accreting neutron stars. In this paper we calculate fusion 
cross sections using Time Dependent Hartree Fock sim- 
ulations coupled with density constrained Hartree Fock 
calculations to deduce an effective potential. We find re- 
markable agreement with experimental cross sections for 
the fusion of stable nuclei. Note that our simulations use 
the SLy4 Skyrme force that has been previously fit to the 
properties of stable nuclei. No parameters have been fit 
to fusion data. 

We compare our results to the simple Sao Paulo barrier 
penetration model. This model calculates an effective po- 
tential by folding over static densities for the projectile 
and target. Within, very roughly, an order of magnitude, 
our results agree with the Sao Paulo model. This pro- 
vides an error estimate for Astrophysical applications of 
the Sao Paulo model. 

However in more detail, there are very interesting dif- 
ferences between our calculations and the Sao Paulo 
model. These differences are likely due to additional dy- 
namics, that is included in our calculations but that is 
neglected in the Sao Paulo model. For the asymmet- 
ric systems ^•^C+^^'O or ^^0-1- ■^''0 we predict an order 
of magnitude larger cross section than for the Sao Paulo 
model. This is likely due to the dynamical effects that 
change the nuclear densities during the collision process 
caused by the rearrangement of the single-particle wave- 
functions. This enhancement of fusion cross sections of 
very neutron rich nuclei can be tested in the laboratory 
with radioactive beams. 

Finally, fusion in accreting neutron stars may take 
place in the inner crust where there is a background gas of 
neutrons. This background gas could impact fusion cross 
sections. To investigate this we are performing TDHF 
simulations where the initial conditions involve two nu- 
clei, appropriately boosted towards each other, and some 
unbound neutrons, that are uniformly distributed in the 
simulation volume. We will report these results in a later 
publication. 
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